function [spectra, descriptions] = ...
  generatePoseFrequency( pose, headers, features, Fs )
%generatePoseFrequency - Computes the pose frequency.

  if ~exist('Fs', 'var')
    Fs = 25;
  end
  
  % Initialize some variables
  N = size(pose, 1);
  spectra = zeros(length(features), N);
  descriptions = cell(length(features),1);
  
  % Iterate through each feature, transforming each into
  % a frequency spectrum
  addpath utils
  for f=1:length(features)
    [Y, desc] = generatePoseFeatures(pose, N, features(f), headers);
    spectra(f,:) = fourierTransform(Y', N, Fs, false);
    descriptions{f} = desc;
  end
  
end

